Sharpening the predictions of big-bang nucleosynthesis 
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Motivated by the recent measurement of the primeval abundance of deuterium, we re-examine 
the nuclear inputs to big-bang nucleosynthesis (BBN). Using Monte-Carlo realization of the nuclear 
cross-section data to directly estimate the theoretical uncertainties for the yields of D, 3 He and 7 Li, 
we show that previous estimates were a factor of 2 too large. We sharpen the BBN determination 
of the baryon density based upon deuterium, ps = (3.6 ± 0.4) x 10~ 31 gem -3 (Q.Bh 2 = 0.019 ± 
0.0024), which leads to a predicted 4 He abundance, Yp = 0.246±0.0014 and a stringent limit to the 
equivalent number of light neutrino species: N„ < 3.20 (all at 95% cl). The predicted 7 Li abundance, 
( 7 Li/H) P = 3.5± ;g x 10~ 10 , is higher than that observed in pop II stars, 1.7 ± 0.3 x 10~ 10 (both, 
95% cl). We identify key reactions and the energies where further work is needed. 
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I. MOTIVATION 

Big-bang nucleosynthesis is an observational corner- 
stone of the hot big-bang cosmology. For more than two 
decades the predicted abundances of the light elements 
D, 3 He, 4 He and 7 Li have been used to test the con- 
sistency of the hot big-bang model at very early times 
(t ~ 0.01 - 200 sec) |,§. The state of affairs in 1995 
was summarized by a concordance interval for the baryon 
density, fl B h 2 = 0.007 - 0.024, for which the predicted 
abundances for all four light elements were consistent 
with the observational data p|. In addition to testing 
the standard cosmology, BBN also gave the best deter- 
mination of the baryon density and was the linchpin in 
the case for nonbaryonic dark matter. 

The big-bang abundance of deuterium is most sensitive 
to the baryon density ||, making it the "baryometer." 
However, deuterium is fragile and is destroyed by stars 
even before they reach the main sequence. Thus, local 
measurements of its abundance, where about 50% of the 
material has been through stars, do not directly reflect its 
primeval abundance. Recently, the situation has changed 
dramatically, ushering in a new, precision era for BBN . 
Buries and Tytler measured the deuterium abundance in 
high-redshift (z > 3) hydrogen clouds where almost none 
of the material has been processed through stars, and 
they have made a strong case for a primeval deuterium 
abundance, (D/H) P = (3.4 ± 0.25) x 10~ 5 [|§. 

From this 10% measurement of (D/H)p, the baryon 
density can be inferred to around 10%, Qsh 2 — 0.019 ± 
0.002, or in terms of baryon-to-photon ratio, rj = (5.1 ± 
0.5) x 10~ 10 . With the baryon density in hand, one can 
predict the abundances of the other three light elements. 
Then, 4 He and 7 Li can test the consistency of BBN; D 
and 3 He can probe stellar processing since BBN; and 7 Li 
can test stellar models. Furthermore, a precise deter- 
mination of the baryon density can make BBN an even 
sharper probe of particle physics (e.g., the limit to the 
number of light particle species) . 

To take full advantage of BBN in the precision era 
requires accurate predictions. The uncertainty in the 
deuterium-inferred baryon density comes in almost equal 
parts from the (D/H) measurement and theoretical er- 
ror in predicting the deuterium abundance. The BBN 
yields depend upon the neutron lifetime and eleven nu- 
clear cross sections (see Table ||). In 1993, Smith, Kawano 
and Malaney (SKM) estimated the theoretical uncertain- 
ties by using a conservative approach ||. While their 
work has set the standard for the past five years, it 
was not without its shortcomings: experiments were not 
weighted strictly by their precision; treatment of system- 
atic effects was neither uniform nor explicit; fits used 
theoretical rules of thumb. In addition, there have been 
new measurements 

After a careful analysis and updating of the micro- 



TABLE I. For each reaction and nuclide, the energies (in 
keV, center-of-mass) at which the sensitivity functions for D 
and 7 Li attain half their maximum value; these intervals in- 
dicate the energies relevant for BBN {Q. B h 2 = 0.019). 



Reaction 


D 


7 Li 


p(n, 7 )d 


25-200 


17-153 


d(p, 7 ) 3 He 


53-252 


65-270 


d(d,p) 3 H 


55-242 


134-348 


d(d,n) 3 He 


62-258 


79-282 


3 He(a, 7 ) 3 Be 
3 He(d,p) 4 He 


no effect 


157-376 


187-325 


107-283 


3 He(n,p) 3 H 
7 Li(p,a) 4 He 


52-228 


24-188 


no effect 


57-208 


7 Li(p,n) 7 Be 


no effect 


1649-1690 


3 H(a, 7 ) 7 Li 
3 H(d,n) 4 He 


no effect 


62-162 


176-338 


167-285 



physics for small but important effects, the theoretical 
uncertainty in the predicted 4 He abundance has been re- 
duced to that in the neutron lifetime, AYp = ±0.001 
(95% cl) Q . Motivated by the primeval deuterium mea- 
surement, we decided to refine the error estimates for the 
other light elements, using the nuclear data themselves 
and Monte-Carlo realization to make our error estimates. 
This method also allowed us to identify where improve- 
ments in the nuclear data would be most useful. 



II. METHOD AND RESULTS 



The details of our method are described in a longer pa- 



per 1 11 1; here we outline the salient points. The nuclear 
inputs come in the form of measurements of cross sec- 
tions, cr{E), or equivalently, the astrophysical S'-factor, 
S(E) = Ea(E)e 2n< > , where e _27r ^ is the Coulomb-barrier 
tunneling probability. From these, the needed thermally- 
averaged reaction rates per particle follow 



(av) 



7T/Lt(fcT) 3 



J a(E)Ee~ E/kT dE, 



(1) 



where fi is the reduced mass. 

We use Monte-Carlo realizations of all the experimen- 
tal data sets to determine thermal reaction rates and final 
yields. For each realization, we proceed as follows. For 
every data point from every data set we draw a value from 
a Gaussian distribution whose mean is the central value 
and whose variance is the standard error reported for that 
point. We account for correlated normalization error in a 
data set by similarly drawing a value for the overall nor- 
malization. For each reaction, a smooth representation of 
S(E) is obtained by fitting a piecewise spline to the data, 
with individual points weighted by their standard errors 
in the usual way. Using the spline fits, we evolve light- 
element abundances with a standard BBN code. From 
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FIG. 1. Summary of the 95% confidence intervals for the 
BBN predictions for D, 3 He, 4 He and 7 Li. The 4 He uncer- 
tainty comes from Ref. Jl^]. Boxes indicate 95% cl abun- 
dances from observation. The vertical band indicates the 
baryon density we infer from the deuterium observation. 



25,000 such realizations, we produce distributions of the 
light-element yields and compute means and 95% cl in- 
tervals. Our results, as a function of the baryon density, 
are shown in Fig. |l|. 

Data points and uncertainties were extracted from a 
comprehensive review of the experimental literature from 
approximately 1945 onward, beginning with a careful 
reading of the original sources. We excluded a small num- 
ber of data sets for which insufficient information for our 
technique was provided. 

As always, there is the sticky problem of systematic 
error, especially for cross sections represented by only 
a few measurements. A case in point is the reaction 
3 He(a,7) 7 Be, which produces nearly all of the 7 Li for 
r] = 5.1 x 10~ 10 . Activation measurements [^2|~H show 
an apparent disagreement with prompt-photon measure- 
ments (see Fig. || and Ref. jlU). Because these mea- 
surements are not in the energy range of relevance for 
BBN, they have little influence on our results. (SKM 
omitted activation measurements from their analysis al- 
together.) We take them into account by performing a 
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FIG. 2. S(E) vs. E C m for the reaction 3 He(a, 7) 7 Be. The 
data are shown with our best fit and 95% cl interval for S(E) 
(solid lines). The SKM fit and 95% cl interval are shown 
as dash-dotted lines. The integration limits on the thermal 
averages needed for an accuracy of O.lcr (inner tick marks) 
and one part in 10 5 (outer tick marks) in the hnal yields are 
shown. g(E) quantifies the sensitivity of the final abundance 
of 7 Li to S(E); the final abundance of D is insensitive to this 
reaction. 

second Monte Carlo, where the prompt-photon measure- 
ments are renormalized by the weighted mean (and un- 
certainty) of the three activation measurements. This 
shifts the 7 Li/H 95% cl interval upward by 11% (see 
Fig. |). 

Our method breaks down completely for the process 
p + n — » d + 7. This is because of a near-complete lack 
of data at the energies relevant for BBN. The approach 
used for this reaction is a constrained theoretical model 
that is normalized to high-precision thermal neutron cap- 
ture cross-section measurements. In particular, we use 
the most recent evaluation, from ENDF-B/VI p6[ . This 
evaluation was performed around 1970 (with a minor up- 
date in 1989), and it fitted a capture model to data of 
similar vintage for the neutron-proton system. No doc- 
umentation survives, and the uncertainty is difficult to 
quantify — especially in light of known systematic prob- 
lems with the likely input data |l7| . (Efforts are under- 
way to construct a new model for this reaction, based 
upon more modern nuclear models and data [[18[ .) For 
consistency, we follow SKM and assign a 5% ler uncer- 
tainty in the overall normalization (also consistent with 
an estimate from the evaluation's authors [0), and we 
use this value for our Monte-Carlo calculations. 

To investigate the role of each reaction independently, 
we ran the BBN code using the SKM rates for all but one 
reaction, studying that reaction alone with our Monte- 
Carlo method. This produced, for each of the eleven 
key reactions, a best fit to the cross-section data, 95% 
cl intervals for the cross sections (Fig. ||), and 95% cl 
uncertainties for D and 7 Li yields for each reaction (see 
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FIG. 3. Uncertainties in the predicted deuterium abun- 
dance from SKM, our full Monte Carlo, and individual re- 
actions, compared with the Buries & Tytler ffl measurement. 
The uncertainties due to reactions not shown are not impor- 
tant. 



ri = 5.1 x 10 



} I i 



& 4 4 



** E ^ 
^ ^ ^ <^ 



FIG. 4. Same as Fig. |, but for 7 Li. The results for the 
alternative normalization of 3 He(a,7) 7 Be are also shown. 



Figs. | and @). 

Our most important result is apparent: the uncertainty 
estimate from our method is a factor of two smaller than 
the SKM estimate. Not only have we reduced the theo- 
retical error estimate, but we have also put it on a firmer 
footing. Our "most probable" yields also differ slightly 
(less than la) from the corresponding results of SKM. 
This reflects both differences in weighting the nuclear 
data, and the inclusion of new data. 

We computed "sensitivity functions" for the yields of 
D and 7 Li for each reaction. These functions measure 
the fractional changes in yield caused by a delta-function 
change in cross section at a given reaction energy (see 
Fig. ^ and Table |). The sensitivity functions quantify 
where precise cross section measurements are required. 



III. DISCUSSION AND CONCLUSIONS 

We have reduced the theoretical error estimate for 
BBN deuterium production by a factor of two, so that the 
deuterium abundance itself dominates the uncertainty 
in baryon density. The deuterium determination of the 
baryon density is thus sharpened, from 8% to 6% (at la), 
or n B h 2 = 0.019±0.0024 (95% cl). In the next five years, 
the precision of the primeval deuterium measurement 
should improve significantly, because the Sloan Digital 
Sky Survey will increase the number of QSOs with mea- 
sured redshifts by a factor of almost 100, with a similar 
increase in the number of deuterium systems expected. 
Further improvement in the theoretical prediction is pos- 
sible; the key reactions in this regard are: d(p,7) 3 He; 
d(d,p) 3 H above 100 keV; d(d,n) 3 He above 100 keV; and 
p(n, 7 )d at 30-130 keV (see Fig. |). Turning the deu- 
terium determination of the baryon density into a few 
percent measurement will make possible a beautiful con- 
sistency test || : comparison with a similarly accurate 
measurement of the baryon density from CBR anisotropy. 

The deuterium-inferred baryon density leads to a pre- 
diction for the big-bang 4 He abundance: Yp = 0.246 ± 
0.001 (D/H) ± 0.001 (t„) = 0.246 ± 0.0014 (all 95% cl). 
When the primeval 4 He abundance is determined to three 
significant figures, this will be a powerful consistency test. 
At the moment, systematic effects dominate the error 
budget; in particular, underlying stellar absorption in the 
most metal-poor Hn regions. Izotov and Thuan's sam- 
ple JT9| , po[ excludes the tainted or suspected-to-be tainted 
systems, and they find Yp — 0.244±0.002. This is consis- 
tent with the deuterium prediction. A less homogeneous 
sample pjj , which includes some of the tainted systems, 
indicates a lower value, Yp = 0.234 ± 0.002, which is not 
consistent with the deuterium prediction. 

Additional light particle species present around the 
time of BBN lead to increased 4 He production, and an 
upper limit to the primeval 4 He abundance can be used to 
constrain their existence p2j. Using Y P = 0.244 ± 0.002, 
the deuterium-determined baryon density, and the prior 
N v > 3.0, we derive the 95% cl limit, N v < 3.20. One 
should be mindful that systematic error in Yp could 
change the limit, and that it will become more secure 
with better 4 He measurements. 

Finally, we turn to 7 Li, the light clement for which 
the uncertainty in the predicted abundance is largest. 
Our analysis has reduced the theoretical uncertainty by 
a factor of two, though a small systematic uncertainty 
remains. Using our full Monte Carlo with the deu- 
terium observations, we predict ( 7 Li/H)p = [3.5±J;g + 
0.4 (sys)] x 10~ 10 . The abundance derived from old, 
pop II halo stars is ( 7 Li/H) pop n = [1-73 ± 0.1 (stat) ± 
0.2 (sys)] x 10- 10 = (1.73 ±0.3) x lO" 10 || (all at 95% 
cl). The discrepancy could represent a real inconsistency 
or merely a depletion of 7 Li by a factor of around two 
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in these stars (predicted by some models of stellar evo- 
lution [QUI). A nuclear solution for the discrepancy 
is unlikely — a ~ 25% (or 5er) change in the p(n, 7)D 
rate would be required, and the unresolved systematics 
of 3 He(ct, 7) 7 Be can only make the problem worse. There 
is still much room to improve the BBN 7 Li prediction; 
the key reactions are: p(n,j)d; 3 He(a,7) 7 Be; d(d, n) 3 Hc; 
and d(p, 7) 3 He. 

Perhaps the most rewarding result of this work is that 
we have verified what David Schramm many times pro- 
claimed, "the predictions of BBN are very robust because 
the key cross section are measured at the energies where 
they are needed." In particular, if all eleven critical cross 
sections were set to zero outside the intervals where they 
are measured, the final light-element abundances would 
change by less than 10% of their current theoretical un- 
certainty. 
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